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We present an analytic five-loop calculation for the additive renormaliza- 
tion constant A{u,e) and the associated renormalization-group function B{u) 
of the specific heat of the 0(n) symmetric cp'^ theory within the minimal 
subtraction scheme. We show that this calculation does not require new five- 
loop integrations but can be performed on the basis of the previous five-loop 
calculation of the four-point vertex function combined with an appropriate 
identification of symmetry factors of vacuum diagrams. We also determine 
the amplitude functions of the specific heat in three dimensions for n = 1, 2, 3 
above Tc and for n = 1 below Tc up to five-loop order. Accurate results are 
obtained from Borel resummations of B{u) for n = 1,2,3 and of the amplitude 
functions for n = 1. Previous conjectures regarding the smallness of the re- 
summed higher-order contributions are confirmed. Borel resummed universal 
amplitude ratios A^/A^ and a^/a^ are calculated for n = 1. 
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I. INTRODUCTION 



One of the fundamental achievements of the renormahzation-group (RG) theory of crit- 
ical phenomena is the identification of universality classes in terms of the dimensionality d 
of the system and the number n of components of the order parameter . Specifically, RG 
theory predicts that the critical exponents, certain amplitude ratios and scaling functions 
are universal quantities that do not depend, e.g., on the strength of the interaction or on 
thermodynamic variables (such as the pressure). The superfluid transition of ^He belongs 
to the d = ?>, n = 2 universality class and provides a unique opportunity for an experimen- 
tal test of the universality prediction by means of measurements of the critical behavior at 
various pressures P along the A-line T\{P). Early tests have been performed by Ahlers and 
collaborators and consistency with the universality prediction was found within the experi- 
mental resolution 0. At a significantly higher level of accuracy, the superfluid density and 
the specific heat (or, equivalently, thermal expansion coefficient) above and below Tx{P) are 
planned to be measured in the Superfluid Universality Experiment (SUE) P] under micro- 
gravity conditions or at reduced gravity in the low-gravity simulator ||^. As demonstrated 
recently 0, this would allow to perform measurements up to |t| ~ 10^^ in the reduced 
temperature t = {T - Tx{P)) /Tx{P). 

On the theoretical side, the corresponding challenge is to calculate as accurately as 
possible the properties of the 0{n) symmetric 0^ model in three dimensions. To extract the 
leading critical exponents from the experimental data and to demonstrate their universality 
at a highly quantitative level requires detailed knowledge on the ingredients of a nonlinear 
RG analysis 0. They include not only the well-known RG exponent functions of the 0^ 
model whose fixed point values determine the critical exponents but also the less well-known 
amplitude functions [[^-[T2[| which contain the information about universal ratios of leading 
and subleading amplitudes 

The existing theoretical predictions on the critical exponents within the minimal 



subtraction scheme p^,p!5| are based on field-theoretic calculations to five-loop order [|T^-[T9 
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and Borel resummation. By contrast, the present theoretical knowledge of the amplitude 
ratios for n > 1 below is based only on low-order (mainly 1- and 2-loop) calculations 
which imply an uncertainty at the level of at least 10-30% [|1|. It has therefore been proposed 
pO| to significantly reduce this uncertainty by performing new higher-order field-theoretic 
calculations and Borel resummations of various amplitude functions in three dimensions. 

Both conceptual and computational steps towards this goal have already been performed. 
The conceptual progress includes the demonstration that the d = 3 field theory suggested 
by Parisi |2ll] can well be realized within the minimal subtraction scheme at d = 3 



by incorporating Symanzik's non-vanishing mass shift p2[ and that spurious Goldstone 
singularities for n > 1 below can well be treated within this approach [0 by using an 
appropriately defined pseudo-correlation length 0. The computational steps include the 
determination of the amplitude functions F^{u) and F^{u) of the specific heat in three 
dimensions above Tc for n = 1,2,3 [|T0| and below Tc for n = 1 [|ll|, respectively, up to 



five-loop order, and their Borel resummation. These calculations, however, were not yet 
complete because of an approximation regarding the additive renormalization A{u, e) of the 
specific heat and the associated RG function B{u). Due to the lack of knowledge in the 
literature about higher-order terms, A{u,e) and B{u) were approximated by their two-loop 
expressions. Although the good agreement between low-order d = 3 perturbation results 
||7| j23|j2^ and accurate experiments |0,|25|j2^ provided some indication for the smallness of 



the effect of the higher-order terms of B{u), no reliable estimate could be given for the 
remaining uncertainty of -F±(m) which could well be of relevance at the level of accuracy 
anticipated in future experiments 0. Furthermore we recall that any inaccuracy of B{u) 
enters not only the formulas ||^ for several universal amplitude ratios but also the formulas 
needed to determine the effective coupling u{l) from the specific heat |7| J23|J2^ . 



It is the purpose of the present paper to provide the missing information on the higher- 
order terms of A{u,6) and B{u) by means of a new five-loop calculation. We shall show 
that the analytic calculation of A{u,e) and B{u) can be directly related to the previous 
calculations [[TBHTSl] of the four-point vertex function. This provides the crucial simplification 



that no new evaluations of three-, four- and five-loop integrals are necessary but that only 
a new determination of symmetry and 0{n) group factors of vacuum diagrams is sufficient. 

Using the five- loop expression of B{u) we are in the position to determine the correct 
higher-order terms of the minimally renormalized amplitude functions for n = 1, 2, 3 

and F_ (u) for n = 1 on the basis of previous work P7|^D| where a different renormalization 
scheme was used. The new coefficients of the higher-order terms of F+{u) turn out to differ 
considerably from the previous approximate coefficients whereas the coefficients of -F-(m) 
are only weakly affected by the new higher-order terms of B{u). 

We also perform new Borel resummations of B{u) for = 1,2,3 as well as of F-{u) 
and of F_{u) — F^{u) for n = 1. It turns out that the result of the Borel resummation for 
B{u) including the new terms up to five- loop order differs from the two- loop result B{u) = 
n/2 + 0{u^) by only less than 1% at the fixed point. For the amplitude functions, our new 
Borel resummation results differ from the previous ones [p!0| , p!I| by about 1% for F_ and by 
less than 0.1% for F_ — -F+ at the fixed point. This is a nontrivial and important confirmation 
of the previous conjectures about the smallness of resummed higher-order contributions 
P|^TT|. As a first application, we calculate the universal ratios A~^/A~ and /a~ of the 
leading and subleading amplitudes of the specific heat above and below T^. for n = 1 and 
c? = 3. Fully quantitative calculations for n > 1 have to be postponed until higher-order 
results for F_{u) are available. 

II. ADDITIVE RENORMALIZATION OF THE SPECIFIC HEAT 

The 0{n) symmetric 0^ model is defined by the usual Landau-Ginzburg- Wilson func- 
tional 

H{U^)} = d'x Qro0^ + i E(V0oi)' + uo{<Piy - ho ■ 0o) (2.1) 

for the n-component field 0o(x) = (0oi(x), . . . , 0on(x)) where 

r-o = r-oc + aot, t = {T - T,)/T, , (2.2) 
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and Hq = (ho, 0, . . . , 0). The Gibbs free energy per unit volume (divided by fc^T) is 

Fo(ro, uo, ho) = -V'^ In J P0o expi-H) . (2.3) 

We shall consider the bulk limit V ^ oo. We are interested in the specific heat per unit 
volume at vanishing external field ho = (divided by Boltzmann's constant ks) where ± 
refers to T > Tc and T <Tc, respectively. Near Tc, (7^ is determined by M 



52 ^ 

„2- 



C^ = Cb- T2^Fo(ro, uo, 0) = Cb - ao^i^o(ro, Uo, 0) (2.4) 



where is an analytic background term. Alternatively the Helmholtz free energy per unit 
volume ro(ro, uo, Mq) = F{ro, uo, ho) + hoMo with Mq = (0oi) determines in the ho ^ 
limit according to 

d^ 

= Cb- al—^To (ro, Mq, Mo(ro, %)) • (2.5) 
org 

The perturbative expression for ro(ro, uo, Mq) is obtained from the negative sum of all one- 
particle irreducible (1 PI) vacuum diagrams |jl5|. The perturbative expression for is 



then determined by the vertex functions f|^'°'* = d^Fo/drQ which we consider as functions 
of appropriately defined correlation lengths ^+ and ^_ above and below Tc PJT^, 



C^ = CB-alt^2''\^i,Uo,d). (2.6) 

A description of the critical behavior requires to turn to the renormalized vertex functions 

ri''°) (e±, u, /i, d) = Z^f (e±, ix'Z^Z^'A,\ d) - ^ix-'A,Aiu, e) . (2.7) 

We work at infinite cutoff using the prescriptions of dimensional regularization and minimal 
subtraction at fixed dimension 2 < d < 4 without employing the e = 4 — d expansion f^-Q . 
The Z-factors are introduced as 

r = Z;\ro-roc), u = ^I-'A,Z-' Z^uo, = V^Vo (2.8) 
where the geometric factor 



Ad = r(3 - d/2) 2^-'^7c-'^/\d - 2)-i 



(2.9) 



becomes = (47r) ^ for (i = 3 and A4 = (Svr^) ^ for c? = 4. These Z-factors Zi{u,e) and 
the associated field-theoretic functions [HI 



C^{u) = /i (9^1nZ<^(u,e)~^ ^ , 
Pu{u, e) = -ew + /3{u) = u\-e + fi (9^(Z^^Zj' 

are known up to five-loop order [[r^ -[l9(|. 



(2.10) 
(2.11) 
(2.12) 



The main quantity of interest in the present paper is the renormalization constant A{u, e) 
in Eq. ( p.TD which absorbs the additive poles of both f and f Previously |TO|JTT[| 
A{u, e) was employed only in its two-loop form 



u 



A{u, e) = -2n- - 8n{n + 2)— + 0{u 



(2.13) 



e 

Here we report on a calculation of A{u, e) up to five- loop order. We would like to stress that 
this calculation does not require new five-loop integrations but can be performed on the 
basis of the previous five- loop calculation [17-0 of the four-point vertex function combined 
with an appropriate identification of symmetry and 0(n) group factors of vacuum diagrams 
which are shown in Fig. |l]. Their negative sum determines the Helmholtz free energy Fq 
up to five-loop order. For the present purpose of determining the pole terms at = 4 it 
suffices to consider only the case tq > and Mq = where only four-point vertices exist. 
The diagrams are labeled (1) in one-loop order, (2) in two-loop order, (3), (4) in three-loop 
order, (5)-(8) in four-loop order and (9)-(18) in five-loop order. The analytic expression 
of an m-loop diagram (i) is given by the product of the coupling (— mo)™'~^, the symmetry 
factor the O (n) group factor G^'^\n) and the momentum integral expression I^^\rQ,e). 
Thus the structure of the diagrammatic expression of a typical diagram, e.g., (16) is 




(2.14) 



2592 (-uo) 



4 + 8n^ + 32^2 + 40n 



X / / / / / Gi ■ G2 ■ Gi+2-3 ■ G3 ■ Gi+2-4 ■ G4 ■ 6*1+2-5 ■ 



(2.15) 



with /p = (27r)~ / d p and the propagators Gi±j = (ro + |pi ± Pj| )~ . The symmetry and 
group factors are hsted in Table |. 

To calculate the additive renormalization constant A{u, e) one needs to calculate those 
ultraviolet d = 4 pole terms of the diagrams contributing to T^q^^ that are left after sub- 
traction of subdivergences. One obtains Fq^^'' by taking two derivatives of Tq with respect 
to Tq. The analytic calculation of the poles of the diagrams for Fq^^'' is identical to that car- 



ried out previously [|T^-[T^ for the four-point vertex function Fo°'^\ To see this, one should 
take into account that in the minimal subtraction scheme the ultraviolet pole terms 
specified above do not depend on tq. Then by using the method of infrared rearrangement 
01 one can nullify ro and introduce for each diagram a new fictitious external momentum 
to regularize infrared divergences. Then one can see that only a particular subset of those 
diagrams of Fq"'^'' are relevant in the present context, namely those where the four external 
legs are connected to each diagram through only two four-point vertices (rather than three 
four-point vertices or four four-point vertices). The details of the calculation are presented 
in Appendix A. 
The result reads 



m=l 

where A^™'^ denotes the contribution of m-loop order 



u 



(2.16) 



A(3)(m,, 



n{n + 2) 



3 40 ^ 24(n + 4) 



u 



(2.17) 



A^^\u,e) = --n{n + 2) 
3 



(n + 8)(12C(3) -25) 96n + 696 



248^ + 1024 48(n + 4)(n + 5) 



n(n + 2) 

15 ^ ' 



u 



(2.18) 



768(n + 4)(n -t- 5)(5n + 28) 128(293^2 + 2624n + 5840) 



9216C(3)(5n + 22) + 32(519^2 + 8462n + 25048) 

192C(3)(7n2 - 28n + 48) + 4608C(4)(5n + 22) + 64(31n2 + 2354n + 9306) 



+ (48C(3)(3n^ - 382n - 1700) + 288C(4)(4n^ + 39n + 146) 

u^, (2.19) 



-3072C(5)(5n + 22) - 3(319^^ - 13968n - 64864)' ^ 



e 



where ({s) = Y.f=i3~' is the Riemann zeta function with C(3) = 1.20205690, C(4) = 7r^/90 
and C(5) = 1.03692776. Most important is the (i-independent RG function B{u) which is 
determined by 

AB{u) = [2Uu) - e] A{u, e) + P^{u, e)^^^ • (2.20) 



Using A{u, e) of Eqs. ( p.l6| )-( p.l9| ) and the perturbative expressions for (r and (3u of Refs. 
iMl we find 



n 8 
B{u) = - + 3n{n + 2)u'^ - -n[n + 2)(n + 8)(25 - 12C(3))m3 
2 3 

+ in(n + 2) [l6C(3)(3r22 - 382n - 1700) - 1024C(5)(5n + 22) 

+96C(4)(4n2 + 39n + 146) - 319^^ + 13968ra + 64864] 

+ 0(m^). (2.21) 

In Table |I| the coefficients c^^ of the power series 

oo 

B{u) = Y: c,^u^ (2.22) 

m=0 

are given for n = 0,1,2,3 up to m = 4 corresponding to five-loop order. Table y also 
contains the coefficients f^"^^ of the power series of the functions 

oo 

m = E (2-23) 

m=l 

where fi{u) denotes the functions P{u), Cr(w) and C(f,{u) for i = 1,2,3, respectively. These 
coefficients are taken from Refs. [p!7|-[T9|l. Up to four-loop order they agree with those in 
Table 1 of Ref. [g. (Note that /f ^ in the table caption of Ref. § should read f^^ x 10"^.) 
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The five-loop coefficients fl^\ /2^'' and /3^'' differ from those in Table 1 of Ref. according 
to the corrections in five-loop order in Ref. [0. 

In Fig. 1^ the different partial sums of B{u) from two- to five- loop order are shown for 
the example n = 2. As expected, the contributions for m > 2 have alternating signs and 
increase considerably in magnitude. Clearly a resummation of B[u) is necessary similar to 
that for Cr{u), C(j,{u) and (3{u) performed previously 

First we reexamine the fixed point values m*, i3u{u*, 1) = 0, for n = 1, 2, 3 obtained in 
Refs. |]8|,|2^ by means of Borel resummation on the basis of previous five-loop results |T7|JT8|] 
and in Ref. [|l^ on the basis of four-loop results. Here we employ the corrected five-loop 
coefficients for the e expansion of the fixed point value which we have derived from Eq. (8) 
of Ref. |jT9| (see Table |ITTD . Employing the standard Borel resummation method we 



have obtained the fixed point values in three dimensions 

u'' = 0.0404 ± 0.0003 forn = 1 , (2.24) 

u* = 0.0362 ± 0.0002 forn = 2 , (2.25) 

= 0.0327 ±0.0001 forn = 3. (2.26) 



The corresponding resummation parameters a and b = 5.5 + n/2 |^ are 

2.22 <a< 3.41 , 6 = 6.0 for n = 1 , (2.27) 
2.45 <a< 3.43 , b = 6.5 for n = 2 , (2.28) 
2.71<a<3.43, 6 = 7.0 for n = 3 . (2.29) 

The previous fixed point values [§,0,^] are consistent with Eqs. ( |2.24| )- (|2.26| ) within the 



previous error bars. The present error bars are smaller than the previous ones [p|, p!6| , p4 
(The range of a determines our error bars, as described further below.) 

We have performed Borel resummations of B{u) at the fixed point u* for the cases n = 1, 
2, 3. In addition, for the important case n = 2 (superfluid ^He), we have determined the 
Borel resummed function B{u) at various values of u. The results are given in Eq. ( [2.341 )- 



A description of the Borel resummation method for the present purpose has been 
given in Section 5 of Ref. (In Eq. (5.10) of Ref. ^ ajk should read aj_rn,k-) In the 
present work, however, we use a different way of determining the parameters a and b of the 
summations. This imphes a different determination of the error bars. 

For B{u) the value of the parameter b is not known from an analysis of the large-order 
behavior (see Eq. (5.6) of Ref. [|] and references therein). Here we fix both b and a by 
requiring fastest convergence of the series of the partial Borel sums S^^^ = S^^\u,a,b) for 
B{u) defined in Eq. (5.12) of Ref. |^ (here L corresponds to (L + l)-loop order). To do so 
we look for the minima of A^^^ and A*^'^^ with regard to variations of both b and a where 



(2.30) 



This yields five-loop values of the parameters a, b for each u. In order to define an error bar, 
we apply the same method to the four- loop result of B{u), i.e. to A*^^) and A*^^\ The four- 
loop values of a, 6 together with the five-loop values provide the ranges of the best values of 
a and 6 as a result of the combined four- and five-loop analysis. Then we define the error 
bar of the five-loop result for B{u) by the maximum and minimum of the resummed four- 
and five-loop values for B{u) over the ranges of the best values of a, b. At the fixed point 
u*, Eqs. (|2.24|) - (|2.26|) , we find the ranges 



0.95<a<1.08, 5.7<6<7.75 forn = 1 , (2.31) 

0.94 <a< 1.04 , 7.0 < 6 < 8.59 for n = 2 , (2.32) 

0.94<a<1.02, 8.18 <6< 9.76 for n = 3 . (2.33) 

The corresponding Borel resummed results for B{u*) are 

B{u*) = 0.5024 ± 0.0011 for n = 1 , (2.34) 

B{u*) = 1.0053 ± 0.0022 forn = 2 , (2.35) 

B{u*) = 1.5080 ± 0.0034 for n = 3 . (2.36) 

We have also determined the function B{u) for n = 2 (superfiuid '^He) in the range < m < 
0.04 as shown in Fig. 
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Most remarkable is the smallness of the deviation of the resummed function B{u) for 
n = 1,2,3 from its two-loop approximation n/2. This confirms previous conjectures [p|-p!l[] 
and justifies earlier analyses 



III. AMPLITUDE FUNCTIONS F± IN THREE DIMENSIONS 



By means of the renormalized vertex functions in Eq. ( p.7|) we define the dimensionless 
amplitude functions 

F±{fi^±,u,d) = -VA^iri''°)(e±,M,/i,rf). (3.1) 

They enter the critical behavior of the specific heat in three dimensions in the form of the 
functions 

F±(1,m,3) = F±(m) (3.2) 

according to 

1 Mi±) 9r (v') — ^ 

C^ = Cb + -a'ii-'A,K^{u{l^)) exp / ^. du' (3.3) 

4 Ju /3u{u', 1) 

where 

K±{u) = F±{u) - A{u,l) (3.4) 

and a = Zr{u, l)^^ao- In Eq- (|3^), u{l) is the effective coupling satisfying 

l^ = PuHl),l) (3.5) 

with u{l) = u. The flow parameters /+ and /_ are chosen as /+ = (/i.^+)^^ and /_ = (/x^_)^^ 
above and below T^. 

The amplitude functions are expandable in integer powers of m and have the power 
series [PHIl 



oo 



m=0 
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and ITT 



Fju) 



u 



E 

m=0 



Cp^U 



(3.7) 



We have determined cj^ up to five-loop order (i.e., up to m = 4) for n = 1, 2, 3 and c~p^ 
up to five- loop order (i.e., up to m = 5) for n = 1 in two different ways. 

(i) The coefficients Cp^ and c~p^ can be calculated from Eqs. (|0|), ( P?7| ) in three dimen- 
sions according to 



(3.8) 



where the Z-factors have the arguments Zi{u, 1). The perturbative expression for F^'^-* can 
be obtained for n = 1, 2, 3 from 



1 



(3.9) 



where the renormalization factor Z^{X) and its relation to the specific heat have been pre- 
sented by Bervillier and Godreche |2^ and by Bagnuls and Bervillier p9| , |32| , see also Ref. 
|T0| . For (i = 3 their renormalized coupling A is related to our uq via the renormalization 
factor Z^{\) according to Mo^+ = — 27rAZ3(A)^^/^ as noted in Ref. [jlO[- The perturbative 
expression of f for n = 1 can be determined according to 



fL''°)(e-,«o,3) 



dr'n 



F_(^_,uo,3) 
1 



d 



d 



F„(^_,Mo,3) 



(3.10) 



where rQ(^_,Mo) is given by Eq. (3.8) of Ref. |ll|]. The Helmholtz free energy f_(^_,uo,3) 
is given by Eq. (3.15) of Ref. |jTl| where our f_(,^_,Mo,3) is denoted by f_o(^-,Mo)- Our 
numerical results for c^^ up to nine digits are presented in Table 

(ii) Alternatively the coefficients can be determined via the relation 

8A^'P±{u)ft'\u) = (1 - 2Cr{u))F±{u) + AB{u) - f3MdF±{u)/du 



(3.11) 



as done previously []TU|,|rT|. For the definition of P± and f±'^^ and for a derivation of Eq. 
( p.ll| ) we refer to Refs. [§] and [§. In the present context we need the contributions to P± 
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and f±'^'' only up to 0{u'^) as given in Table 4 of Ref. |10] and Table 3 of Ref. [11| since B{u) 
is known only up to O(m^) as well. (We recall that the coefficients of P_ are determined by 
those of P+ according to P_(m) = — |{1 + 2[1 — P+{u)] — ^(r{u)} 0.) This calculation via 
Eq. yields coefficients c^^ that agree with those obtained via Eq. ( p.8| ) up to eight 



digits for Cp^ and up to seven digits for Cp^. We consider the calculation (i) via Eq. (|3.8|) as 



slightly more reliable since fewer numerical operations are required than in the calculation 



(ii) using Eq. (KW) 



Since here we have used the perturbative contributions of B{u) up to five-loop order, 
the resulting higher-order coefficients given in Table |I^ differ from those determined 



previously (see Table 4 of Ref. and Table 3 of Ref. |TT[|) where the approximation 



B{u) = n/2 + Oiv?) was used. Only our low-order coefficients CpQ, Cp^, CpQ, Cp^ and Cp2 
agree with the previous ones [|l^,|ll|]. The coefficients with m > 1 differ considerably 
from the previous ones whereas the coefficients Cpj^ with m > 2 differ only by 0.2% (m = 3), 
0.1% (m = 4), and 2% (m = 5). 

In order to study the effect of the new higher-order terms we have performed Borel 
resummations of the series for uF_{u) and for u[F_{u) — F+(u)] at the fixed point u* for 
the case n = 1. The method employed is the same as for B[u) in Sect. II. The parameter 
ranges turn out to be 

1.6 < a < 1.7 , 7.48 <h< 8.70 (3.12) 

for u*F_{u*), and 

1.4<a<1.7, 6.0 <6< 11.7 (3.13) 

for u*[F„(m*) -F+(u'^)]. 

We have found that our present method does not yield a reliable estimate of the param- 
eters a and b for F+(m*) separately; this may be related to the fact that, unlike Cp^^, the 
coefficients (see Table |rV|) do not have alternating signs for m < 3 (this alternation 
is predicted for the asymptotic large-order behavior |]8|,|3l|). P+('u) will be further studied 
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elsewhere. In the apphcation to amphtude ratios given below we shall not need F+{u*) 
separately. 

The resummation results are 



(m*) = 0.3687 ± 0.0040 (3. 14) 

and 

u''[F^{u*) - F+{u*)] = 0.4170 ± 0.0036 . (3.15) 



The previous approximate resummation results |[TO| , p!Tl| forn = 1 were u*F^{u*) = 0.3648 and 



u*[F_{u*) —F^{u*')] = 0.4170 with an error bar of about 1%. Thus our resummation results 
differ from the previous ones only by about 1% for F_ and by less than 0.1% for F„ — 
confirming previous conjectures |10| , p!I[| . The parameter dp in the effective representation 
TT| F^{u) = {2u)~^ — 4(1 + dpu) now becomes dp = —4.64 (compared to —4.04 in Ref. 



Tl|). This leaves the solid line in Fig. 4 of Ref. [|Tl[] essentially unchanged. 

As an illustration we apply our results to the specific heat in the asymptotic critical 
region where it can be represented as 

A± 

C^ = — \t\-'' (l + at\t\^ + ■■■)+ B (3.16) 

a ^ ' 

with the Wegner exponent |3^ A. We consider the universal amplitude ratios [|l| j A~ 
and /a~. The former can be expressed in terms of B* = B{u*) and F^ = F±{u*) in three 
dimensions as iPl 



A- U_ 



F 

1 + a 



AvB* + aF* 



(3.17) 



where 



2uPl 

(3.18) 



6_ (3/2) - 2vPX ■ 

with P| = P+(m*). For the Borel resummed value of P| we have obtained 0.757 ± 0.004 for 
n = 1 . The expression of a+ / a~ is more complicated and depends on the derivatives of the 
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functions F_, F_ — F+, P^, B, and Cr at the fixed point, as specified in Eqs. (4.24), (5.19) 
and (5.21)-(5.24) of Ref. 0. We have performed Borel resummations for these quantities 
on the basis of our new resuhs. Due to the present lack of knowledge on higher-order terms 
of -F_ below Tc for n > 1 we have confined ourselves to the case n = 1. 

Our result for A'^/A^ depends on the value employed for the critical exponent u = 
(2 — a)/3. For u = 0.6310 ||l3l we obtain A'^/A'^ = 0.539, with an uncertainty of about 
2%. For a^/a~ we obtain the value 1.0, with an uncertainty of 0(10%). A more detailed 
presentation of these applications including error bars and a comparison with previous results 
will be given elsewhere. Here we only note that B* and enter also various other important 
universal ratios, e.g., Rj and a~ /ap^ related to the superfluid density ps 0. A quantitative 
calculation of these ratios must be postponed until higher-order results are available for 
n > 1 below T^. 

After completion of the present work we learned of the preprint [hep-ph / 9 71 0346| by B. 
Kastening, "Five-Loop Vacuum Energy Function in 0^ Theory with 0(n)-Symmetric and 
Cubic Interactions" where the perturbative terms of a function equivalent to B{u) have been 
calculated up to five-loop order. These terms agree with ours in Eq. (|2.21| ). 
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APPENDIX A: DERIVATION OF THE ADDITIVE RENORMALIZATION A{u, e) 



The following expressions 



give the ultraviolet d = 4 poles of the diagrams 



shown in Fig. ^defined by K R' {d'^ I^'^^ / dr^) according to the standard notations, see e.g. Ref. 
|TE|. Here R' denotes the incomplete ultraviolet i?-operation which subtracts sub divergences, 
and K denotes the operation of taking pole parts. We use subscripts (a.b) for the pole terms 
on the right-hand-sides of (|AT|)- ([A8|) that are identical with the numbers associated with the 
three- and four-loop diagrams of Ref. |jT6| (the first number a in the brackets indicates the 
number of loops and the second number b indicates the consecutive number of a diagram in 
Table 1 of Ref. [|16|). The subscripts on the right-hand sides of Eqs. ( |A9|) -( [AT^ ) correspond 



to the numbers of the five- loop diagrams of Ref. . We have multiplied the left-hand sides 
by factors (IGvr^)"* in m-loop order because of the definition of the bare four-point coupling 
Wn'^ go/24: in Refs. 

The one- and two-loop expressions read 

1 



(o2 



-- J ln(ro 



{IGnyKR' 
the three-loop expressions read 

/ 02 



drn 



21 



(2.2) 



{16nr KR' /(^)(ro,6) 
{16nrKR' /(^)(ro,6) 



21, 



(3.2) 



J(3) 



8/(3.4) + 12/(3.9) = J^'^ 



the four-loop expressions read 

/ 92 



(4.5) = J^'^ 



aG^rKR'l^^^I^'\ro,e)j =2 

(16vr^r/r/?'(|^/(^)(ro,e)) =0 

{WtcYKR' (^I^ /(''(ro, e)^ = 4/(4.9) + 6/(4.11) = J^'^ 

{lQ7r')^KR' (^1^ /(8)(ro, e)^ = 24/(4.12) + 6/(4.18) + 12/, 



(4.19) 



J(8) 



(Al) 
(A2) 



(A3) 
(A4) 



(A5) 
(A6) 
(A7) 
(AS) 



16 



and the five-loop expressions read 

/ Q2 



/(9)( 



il6.r KR' 

il6.r KR' /(-)(ro,e 
il6nr KR' /(-)(ro,e; 
ilGny KR' (^I^''\ro,e: 



(16.^)^i^i?'(^/(-)(ro,e) 



(167r2firi?'(^/(i«)(ro,e) 



(9ro 
(9ro 



2/nr = J^'^) , 
= , 
= , 
2/121 = , 
4/7 + 6/120 = J^^^^ 
2/8 = , 



(A9) 
(AlO) 
(All) 
(A12) 
(A13) 
(A14) 



2/14 + 12/15 + 12/19 + 24/21 + 4/23 + I8/106 = J^''^ , (A15) 
8/79 + 16/88 + 32/95 + I6/100 = J^''^ , (A16) 

2/81 + 8/84 + 4/99 = /^^), (A17) 

8/32 + 32/47 + 4/73 + 8/93 + 8/98 + 4/109 + 8/116 = J^''^ . 

(A18) 



The pole terms up to four loops are ||T6| 

/ 



'(1.1) 



'(2.2) 



'(3.2) 



2 

e 



2_ _ 3_ 

3e2 ~ 4^ 



/(3.4) = — - — 



'(3.9) 



'(4.5) 



'(4.9) 



+ 



3e3 3e2 3e ' 
16 

4 3 

~3^^ 2^' 
16 16 4 



/(4.11) = ——r + -r^ 



3e4 + 3e3 3e^ 



(A19) 
(A20) 
(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 



17 



'(4.12) 



'(4.18) 



4 10 

3?^ 3? 



13 11-6C(3) 
3e2 6e 



8 8 4 2C(3) 

\ \ h 

3e4 ^ 3e3 e 

1 7 



2 

3^3 



-^(4.19) = -— + — - 



12e 



and the five-loop pole terms are [18 

h 
h 

/l4 
/l5 
/l9 
/2I 
I23 
-^32 
-^47 
-^73 
-^79 
/81 



8 


3 




3f4 ~ 
8 


3 




3e4 ~ 






4 


3 


5 


15e3 


10e2 


~ 96e ' 


2 


13 


857 


15e3 


20e2 ^ 


^ 960e ' 


8 


7 


25 215 


15e4 


■3? + 


6e2 96e ' 


16 


7 


11 157 
60e2 ' 320e 


15e4 


' 5^ " 


4 


3 


5 


15e3 


10e2 


'We' 


48C(3) 


48C(3) + 24C(4) 


5e3 




5e2 


8 


12 

■5? + 


6 18C(3) - 


15e5 


e3 ^ 5e2 



+ 



5e 



+ 



30e 



16 8 
TSe^ ^ 3^ 
16 16 



28 



6 + 4C(3) 32 - 12C(3) - 18C(4) 



5e2 



30e 



8 4 + 4C(3) 7 + 6C(3)-12C(4) 

5e5 5e4 5e^ 5e^ 5e 



16 16 

3? ~ 3? 
20 

3e5 3^ 
16 



8 4 - 4C(3) 



^84 = — - — + — 
-^88 



3e3 e2 
26 44 - 24C(3) 

28 



12e2 

6 - 12C(3) 16 - 18C(4) 



-^93 = 
-^95 

-^98 



15e5 3e4 
52 

5e5 ~ " 
16 



15e4 



+ 



15e3 5e2 15e 

34 116 -24C(3) 56 - 44C(3) + 6C(4) 

15^ ^ 60? We 

28 6 - 12C(3) 16 - 18C(4) 



15e5 3e4 15e3 

4 14 19 

+ 



15e4 15e3 15e2 
4 2 7 

-'99 — ~ T3 + ' 



5e2 15e 
386 + 768C(3) 

960e ' 



3e4 



6e2 



6 



- 5^ " 5^ + 5^ + 



81 -48C(3) 
iWe 



(A27) 
(A28) 
(A29) 

(A30) 
(A31) 
(A32) 
(A33) 
(A34) 
(A35) 
(A36) 
(A37) 
(A38) 
(A39) 
(A40) 
(A41) 
(A42) 
(A43) 
(A44) 
(A45) 
(A46) 
(A47) 
(A48) 



-^106 
Il2l 



64 _ 32 
15e5 ~ 5e4 
16 16 

5^ 
4 

3^ 



15e5 



+ 



15e4 
32 

32 



32 



+ 



5e3 
16 

5^ 
32 

15^ 



3e3 



16 2 

15^ ~ 15e ' 

52- 108C(3) 202 - 168C(3) - 18C(4) 
15^2 

250 - 96C(3) 

120e ' 



3e5 3e4 
32 32 8 

3?~ 3?^3^ 



(A49) 
(A50) 
(A51) 
(A52) 
(A53) 
(A54) 



In Eqs. ( |A37|) and ( |A5 1| ) the corrections found in Ref. ||T9[ have been taken into account. 
Note that in Eqs. ([A19| )-( |A54|) we have used e = 4 — whereas in Refs. |[T6| . |18| e denotes 
(4-rf)/2. 

Eqs. ([A1|)-( |A5^ ) determine the additive renormahzation according to Eq. ( |2.7| ) as 



A{u,t) 



n 8n(n + 2) 

- + ^-2 -i-^ 

6 



i=3 



i=5 



18 
i=9 



(A55) 



The overall factor of 2 in Eq. ( [A55| ) arises from the d = 4 value of the factor (Ai^/A) ^/16 



TT 



which is needed to obtain A{u,e) from Eqs. ( |A1|) -( |A5^) according to Eq. ( p.TI) . The renor- 
malized coupling u enters Eq. (|A55|) in the form u/2; the factor 1/2 takes into account that, 
near d = A, Uq = A^^u + 0(^2) = grr^u + 0{u'^) ^ Wn^u + 0{u^) [see Eq. (^)]. 
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APPENDIX B: Z FACTORS 



In deriving the coefficients of the perturbation series of the quantities -F±('u), P±('u), 
and /±'°''('u) we needed the Z-factors Z.^, Z^, and calculated previously [p!7|-[T9[| up to 
five-loop order. Since their explicit form is not available in the literature we present them 
here explicitly. They read 

5 



Z,[u, e) = 1 + 5] Zf^\e)u'^ + O(m^) 

m=l 
m=l 



(Bl) 
(B2) 
(B3) 



m=l 



with the following coefficients in m-loop order: 
Coefficients of ZJu, e): 





4(n + 2) 
e 


Z?\e) 


= 4(n + 2) 




16 , 
= y (^ + 2) 




2 , 

= — (n + 2 

3 ^ 



4(n + 5) 



15n + 111 



-278 - 61n Vlr? + 132n + 360' 



(B4) 
(B5) 

(B6) 



288C(4)(5n + 22) + 48C(3)(3n2 + lOn + 68) + 31060 - + 7578n 



1152C(3)(22 + 5n) + 1236^2 + 23580n + 74616 ^ 16(245^2 + 2498n + 6284) 



-192 



(n + 5)(2n + 13)(n + 6) 



(B7) 



Zf\^) 



15 



:^+2) 



(9600C(6)(55n + + 186) + 768C(5)(-5n2 + 72 + 14n) 
+ 288C(4)(29n2 + 2668 - 3^^ + 816n) + 768C(3)2(-2n2 - 145n - 582) 
+ 48C(3)(8208n + 17n^ + 940n2 + 31848) + 21n^ + 4525471^ + 1077120n 
+ 3166528) ^ - (30720C(5) {2n^ + 186 + 55n) + 576C(4) (5n + 22) {n - 22) 
+ 96C(3)(27n3 + Vll^n^ + 14456n + 45448) - 9871^ + 277280^^ + 3073376n 
+ 7449712)^ + (2304C(3)(13n + 74) (5n + 22) + 21576^^ + 685192n2 



20 



^o. ^ -. ^ . .r.r..r.\ 1 32(307976 + 3175271^ + 172176n + 1933n=^ 
+ 5017312n + 10459360)— ^ . 



+ 



384(5n + 34) (n + 6) (n + 5) (2n + 13) 



Coefficients of Z^{u, e): 



7(1 
7(2: 
7(3 

7(5 



0, 



4(n + 2) 



= -(n + 2)(n + 8) 



3 

2(n + 2) 



1 4 



5(n2 - 18n - 100) 4(71^ + 234 + 53n) _ 16(n + 8)^ 



-^(» + 2) 



( - 1152C(4)(5n + 22) + 48C(3)(-6n^ + 184 + n^ + 64n) 
- 22752n - 39n^ - 296n^ - 77056)- 



+ 



2304C(3)(5n + 22) - 60n^ + 135488 + 42440n + 1844/1^ 



16(n + 8)(3n2 + 269n + 1210) ^ 192(n + 8)^ 



Coefficients of Zu{u, e): 



^ 16 



(n + SV 5n + 22 



+ 



96C(3) (5n + 22) + 942n + 2992 + 3571^ 16(n + 8) (17n + 76) 



24(n + 8)^ 



zi:^ 



16 

T 



+ 



(480C(5)(2n2 + 55n + 186) - 72C(4)(n + 8)(5n + 22) 

24C(3)(63n^ + 764n + 2332) + 20624n + 1640n^ - 5n^ + 49912)- 
480C(3)(5n + 22) (n + 8) + 67424n + 153088 + 7736/1^ + 172n^ 

16(55n + 248)(n + 8)2 48(n + 8)^ 



4 
15 



(6912C(7) (25774 + 9261n + 686n2) - 28800C(6)(n + 8)(2n2 + 55n + 186) 
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+ 768C(5) (165084 + 7466n^ + 305n^ + 66986n) - 288C(4) (62656 + 4084n^ + 28084n 

+ 189n^) + 2304C(3)^(3264 - 59n^ - 6n^ + 446n) + 48C(3)(1264n^ - 13n^ + 1312864 

+ 551032n + 67432n^) + 20429248n + 2518864n^ + 195n^ + 40148480 + 39230n^) - 

- (99840C(5)(n + S){2n^ + 55n + 186) - 14976C(4)(5n + 22) (n + 8)^ 

+ 3456C(3)(91n^ + 15436n + 34144 + 2196n^) + 63219712n - 800n^ 

+ 420800n^ + 117768192 + 981171271^) 4 

66048C(3)(5n + 22) (n + 8)^ + 32(n + 8)(733n3 + 4018671^ + 353392n + 803328) 

512(193n + 875)(n + 8)3 3840(n + 8)51 ,^ , 
-4 + ■ (B18) 
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FIGURES 

FIG. 1. Vacuum diagrams up to five-loop order determining the Helmholtz free energy Fq for 
Mq = 0, ro > 0. Diagrams (6), (10) and (11) do not contribute to A{u,e). The pole terms derived 
from the vacuum diagrams are given in (|AlD -( A54 ) of App. A. 



FIG. 2. Partial sums Bm{u) = E-m=o Cs^^*"" of B{u), Eq. < ^^ , as a function of u for 



n 



from M = 1 (two- loop order) to M = 4 (five-loop order). Also shown is the Borel resummed 
result (solid line) which deviates from the two-loop result i?i = 1 by only 0.5 % at the fixed point 
u* = 0.0362. 



FIG. 3. Borel resummation result for the function B{u), Eq. ( 2.21| ), for n = 2 (solid line) 



obtained by interpolation between the resummed values of B{u) at = ku* /lO, A; = 1, . . . , 10 of 
the renormalized coupling u in the range < u < = 0.0362, with error bars. Also shown is the 
three-loop result B2{u) = 1 -|- 24n^ (dashed line), compare Fig. |2[ The two- loop result is i?i = 1. 
The Borel values Bk = B{uk) are 1.000233, 1.00073, 1.0013, 1.0019, 1.0026, 1.0032, 1.0037, 1.0043, 
1.0048, 1.0053 for /s = 1, . . . , 10, respectively. 
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TABLES 

TABLE L Symmetry and group factors of the vacuum diagrams shown in Fig. ||. Diagrams 
(6), (10) and (11) do not contribute to A{u,e). 



loop order 


diagram (i) 


symmetry factor 5^*^ 


group factor G^^\n) 


1 loop 


(1) 


1 


n 


2 loops 


(2) 


3 


i(n2 + 2n) 


3 loops 


(3) 


36 


i(n^ + 4n^ + 4n) 




(4) 


12 


i(n2 + 2n) 


4 loops 


(5) 


432 


^{n'^ + 6n^ + I2n^ +8n) 




(7) 


576 


^{n^ + 4n2 + 4n) 




(8) 


288 


^(n^ + lOn^ + 16n) 


5 loops 


(9) 


5184 


^{n^ + 8n'^ + 24n^ + 32n^ + 16n) 




(12) 


10368 


^(n'^ + 6n3 + 12n2+8n) 




(13) 


6912 


^(n'^ + 6n3 + 12n2+8n) 




(14) 


6912 


^(n'^ + 6n3 + 12n2+8n) 




(15) 


2304 


i(n^ + 4n^ + 4n) 




(16) 


2592 


g^(n'' + 8n^ + 32n2 + 40n) 




(17) 


20736 


-^(n^ + 12n3 + 36n2 + 32n) 




(18) 


10368 


g^(5n3 + 32n2+44n) 
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TABLE II. Coefficients /j of the functions f3{u), Cr('u) and C(^(u) for i = 1,2, 3, respectively, 
and coefficients c^^ of B{u), compare Eqs. (2.20) and (2.21), for n = 0, 1, 2, 3 up to five-loop order 



(m = 6 for (3, m = 5 for d^r and C^, m = 4 for B). For compare Table 1 of Ref. 









Cr 




B 


n = 








8 












32 


-80 


-16 









-672 


3552 


128 









43989.9534 


-223152.607 


-8000 









-4166409.19 


18836823.8 


500639.112 









498653403.0 








n = 


1 





12 





1/2 






36 


-120 


-24 









-816 


6048 


216 


9 






56245.8519 


-413813.942 


-14040 


-761.422836 






-5632017.54 


37512804.7 


958294.321 


44244.7100 






708814936.0 








n = 


2 





16 





1 






40 


-160 


-32 









-960 


9024 


320 


24 






69029.7505 


-660870.017 


-21120 


-2256.06766 






-7268274.40 


63662497.1 


1566676.68 


141294.329 






956636505.0 








n = 


3 





20 





3/2 






44 


-200 


-40 






28 



-1104 12480 440 45 

82341.6490 -967074.371 -29000 -4653.13955 

-9075019.76 98265069.9 2333667.84 310944.846 
1243816220.0 

TABLE III. Coefficients fi''^ of the e = A - d expansion of the fixed point value 
u*{e) = Y^k=i fu'^^^ up to = 5 (five-loop order) for n = 1, 2, 3. 

k i!!'' 

n = l 1 1/36 

2 17/972 

3 -0.0114632463 

4 0.0223877393 

5 -0.0703070454 

n = 2 1 1/40 

2 3/200 

3 -0.00896474627 

4 0.0170850033 

5 -0.0492703392 

n = 3 1 1/44 

2 69/5324 

3 -0.00718789532 

4 0.0134615145 

5 -0.0358025667 
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TABLE IV. Coefficients of and for n = 1, 2, 3 defined in Eqs. and 

(|3.7|) , respectively. For cj^„, ?7i refers to corresponding to (m + l)-loop order whereas for Cp^, 
m refers to n™^^ corresponding to m-loop order. 





m 


+ 




n = 1 





-1 


1/2 




1 


-6 


-4 




2 


-22.6976286 


72 




3 


-722.742494 


-5189.75477 




4 


34775.5862 


433582.588 




5 




47754702.4 


n = 2 





-2 


1/2 




1 


-16 


-4 




2 


-92.5270094 


64 




3 


-2430.86459 






4 


102469.659 




77, = 3 





-3 


1/2 




1 


-30 


-4 




2 


-233.488142 


56 




3 


-5742.02974 






4 


204463.778 
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1 loop: 

(1) 



2 loops: 




(5) (6) (7) (8) 



5 loops: 




(9) (10) (11) (12) (13) 




(14) (15) (16) (17) (18) 



B(u) 

1.1 



1.0 



1 and 2 loop 
Borel 


5 loop / 

/ 


n=2 


/ 

/ 

/ 

/ 

/ 

/ 




/ 

Sloop 


1 


4 loop ^ \ 

\ 

1 



0.01 0.02 0.03 0.04 



